model {
	
	
	for (l in 1:L){										
		for (j in 1:N) {
			mrgbar[j,l] <- lambda[l] * phi[mrg.it[j]] 
			mrg.y[j,l] ~ dnorm(mrgbar[j,l], sigma[l])
			}
		}
	
	
	# Phi
	for(j in 1:phi.IT){
		phi[j] ~ dnorm(0,1)
		}
	
	# Lambda / Sigma 
	for(l in 1:L) {
		lambda[l] ~ dnorm(l0, L0)
		sigma[l]  ~ dgamma(e0, f0)
		}	
		
}